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Preface 


This  paper  has  been  written  In  manuscript  format  rather  than  the 
traditional  thesis  format.  The  School  of  Oceanography  encourages  this 
approach  to  reduce  the  work  involved  in  preparing  papers  stemming  from 
graduate  research  projects  for  publication  in  scientific  journals.  For 
this  reason,  some  deviations  from  the  ordering  of  a traditional  thesis 
are  present:  1)  no  list  of  Tables  or  Figures  appears;  2)  Tables  appear 
In  order  at  the  end  of  the  main  text;  3)  a list  of  Figure  legends,  fol- 
lowed by  the  Figures,  are  placed  after  the  Tables;  4)  acknowledgments 
are  placed  after  the  text,  rather  than  before.  This  paper  will  be  sub- 
mitted to  the  Journal  of  Physical  Oceanography  with  Thomas  J.  Spoering 
as  first  author  and  Clayton  A.  Paulson  as  second. 


1.  Introduction 

The  investigation  of  internal  gravity  waves  In  the  open  ocean  has 

produced  an  extensive  theoretical  and  observational  literature.  Recent 

advances  in  the  observations  and  interpretation  have  been  excellently 

reviewed  by  Garrett  and  Munk  (1979).  For  the  most  part,  the  theory  and 

observations  have  applied  to  regions  of  the  ocean  in  which  the  buoyancy 
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frequency,  N,  defined  by  N = g p 3p/Dz,  is  nearly  constant  or  changes 
slowly  with  depth.  In  the  upper  ocean,  N usually  varies  rapidly  with 
depth,  reaching  a maximum  below  the  base  of  the  surface  mixed  layer  and 
then  decreasing  with  depth  in  the  seasonal  thermocline. 

Garrett  and  Munk  (1972,  1975)  successfully  modeled  the  distribution 
of  the  internal  wave  energy  in  wavenumber  and  frequency  space.  This 
model  spectrum  was  based  on  observations  from  many  sources  together  with 
simplifying  assumptions.  The  applicability  of  the  Garrett-Munk  model 
near  the  base  of  the  surface  mixed  layer  is  doubtful. 

Pinkel  (1975)  reported  observations  from  FLIP  between  60  and  400  m 
depth  from  which  spectra  of  isotherm  displacement  and  slope  were  computed. 
Pinkel  found  that  the  first  mode  is  strongly  dominant  above  2 cph  while 
below  2 cph  the  bandwidth  of  energetic  wavenumbers  is  broader,  corres- 
ponding to  multiple  energetic  modes. 

Kase  and  Clarke  (1978)  analyzed  moored  and  profiling  CTD  measurements 
in  the  upper  thermocline  during  the  GARP  Atlantic  and  Tropical  Experiment 
(GATE).  They  found  a peak  in  the  internal  wave  spectrum  in  the  frequency 
band  from  2 to  5 cph.  Waves  in  this  band  are  coherent  and  in  phase  with- 
in the  upper  30  m of  the  thermocline.  Kase  and  Clarke  suggest  a model 
showing  how  surface  wind  stress  forcing  may  produce  a spectral  peak  in 
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the  upper  thermocline  at  frequencies  higher  than  the  buoyancy  frequency 
of  the  thickest  underlying  layer  of  constant  N. 

Nelson  and  Milder  (1979)  analyzed  towed  thermistor  chain  observations 
In  the  upper  150  m of  the  Eastern  Tropical  Pacific.  Spectra  of  isotherm 
depths  agreed  with  the  Garrett-Munk  model  although  there  was  evidence  of 
horizontal  anisotropy  at  wavenumbers  below  1 c/km.  They  found  that  towed 
vertical  coherence  was  substantially  higher  throughout  the  horizontal 
wavenumber  band  when  compared  to  previous  deeper  observations  and  the 
Garrett-Munk  model . 

Dillon  and  Caldwell  (1979)  reported  observations  In  the  upper  200  m 
of  the  vertical  velocity  of  internal  waves  obtained  from  a nearly  free- 
falling  vehicle.  The  spectrum  of  vertical  velocity  has  a maximum  at  a 
band  in  the  spectrum  centered  at  3 cph.  They  found  that  most  of  the 
internal  wave  energy  in  this  band  was  in  the  first  and  second  modes. 

The  objective  of  this  paper  is  to  describe  observations  of  internal 
waves  in  the  upper  North  Pacific.  The  measurements  were  made  by  use  of 
a towed  thermistor  chain  at  depths  between  20  and  40  m.  Horizontal  wave- 
number  spectra  of  the  vertical  displacement  of  isotherms  are  presented 
and  compared  to  previous  observations  and  to  the  Garrett-Munk  model. 
Variability  of  the  spectra  with  the  local  value  of  N,  tow  direction,  and 
wind  speed  are  investigated.  Coherence  spectra  between  isotherms  sepa- 
rated in  the  vertical  are  also  presented. 
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2.  Instrumentation 

The  design  of  the  towed  thermistor  chain  was  based  on  the  design  of 
a similar  system  constructed  by  the  Applied  Physics  Laboratory  of  Johns 
Hopkins  University  (Mobley  et  aK , 1976).  The  main  difference  in  our 
design  is  the  incorporation  of  the  bridge/amplifiers  into  the  chain  near 
each  of  the  thermistors  rather  than  having  the  bridge/ampl ifiers  located 
in  the  ship's  laboratory.  A more  complete  description  of  the  design  and 
construction  of  our  system  is  given  by  Mesecar  and  Evans  (1977). 

A schematic  diagram  of  the  thermistor  chain  in  operation  is  shown 
in  Fig.  1.  The  chain  was  constructed  using  plastic  fairing  manufactured 
by  Fathom  Oceanology  Ltd.  (Model  770T).  The  fairing  effectively  reduces 
drag  and  vibration  and  also  supports  the  sensors,  bridge/amplifiers  and 
electrical  conductors.  A steel  strain  member  runs  through  the  nose  pieces 
of  the  fairing  and  is  attached  to  a streamlined,  lead-filled,  cyclindri- 
cal  depressor  which  has  a mass  of  450  kg.  The  fairing  is  free  to  turn 
around  the  strain  member  and  align  itself  parallel  to  the  flow.  The 
angle  between  the  chain  and  the  vertical  was  less  than  10°  for  tow  speeds 
up  to  3 m/s  with  the  depressor  at  a depth  of  40  m. 

There  were  25  thermistors  installed  on  the  chain,  each  separated  by 
a distance  of  one  meter.  Four  pressure  sensors  and  associated  electronics 
were  installed  at  intervals  of  eight  meters.  The  specifications  for  the 
temperature  and  pressure  measurements  are  given  in  Table  1. 

Sensors  and  associated  electronics  can  be  placed  anywhere  along  the 
length  of  the  chain  in  modified  sections  of  fairing.  The  installation  of 
a thermistor  and  its  associated  bridge/amplifier  is  illustrated  in  Fig.  2. 
Part  of  the  nosepiece  of  a fairing  section  was  cut  out  and  replaced  by  a 


nosepiece  containing  a thermistor.  The  thermistor  was  molded  into  the 
nosepiece  and  exposed  to  the  environment  through  an  indentation  in  the 
front.  The  leads  from  the  thermistor  run  to  another  molding  containing 
pad  resistors  used  to  normalize  the  calibration.  Leads  run  from  the  pad 
resistors  to  a molding  in  the  tail  section  containing  the  bridge/ampli- 
fier. The  bridge/amplifier  is  powered  from  shipboard  by  electrical  con- 
ductors running  through  the  tail  sections  of  the  fairing.  Signals  are 
conducted  to  the  ship  by  similar  means.  The  signal  and  power  conductors 
are  twisted  pairs  of  teflon-coated,  stranded  wire.  The  pressure  sensors 
and  associated  electronics  were  installed  entirely  within  tail  sections 
of  fairing  analogous  to  the  installation  of  the  bridge/amplifiers  in  Fig. 
2. 

A specially  constructed  winch  and  sheave  are  part  of  the  system. 
Because  of  the  stiffness  of  the  chain  and  the  wires  running  through  the 
tail  sections  of  the  fairing,  a large  turning  radius  is  required.  The 
winch  has  a drum  1.5  m in  diameter  which  holds  a single  wrap  of  faired 
cable.  The  drum  has  a continuous  plastic  groove  attached  to  its  surface 
which  is  shaped  to  accept  the  fairing  nosepieces.  The  grooved  surface 
eliminates  the  need  for  a level -wind.  The  chain  runs  out  over  a sheave 
1.2  m in  diameter  which  is  suspended  from  the  ship's  boom  or  A-frame. 

The  signals  from  the  thermistors  and  pressure  sensors  are  processed 
and  recorded  in  the  ship's  lab.  The  signals  are  amplified,  offset, 
digitized  and  recorded  on  magnetic  tape  by  use  of  a minicomputer  system. 
Means  and  standard  deviations  are  computed  simultaneous  to  recording. 

The  thermistors,  together  with  the  bridge/amplifiers,  were  calibrated 
in  a stirred  calibration  bath  which  was  stable  to  better  than  0.001°C. 


Comparison  of  calibrations  before  and  following  the  experiment  showed 
agreement  to  within  ± 0.01 °C  for  most  of  the  sensors.  Disagreement  was 
caused  by  rupture  of  the  glass  insulation  covering  the  thermistor  bead. 

The  calibration  data  were  fit  to  a curve  proposed  by  Steinhart  and 
Hart  (1968): 

T'1  = A + 8 In  R + C (In  R)3 

where  T is  the  temperature  in  degrees  Kelvin,  R is  the  resistance  of  the 
thermistor  and  A,  B and  C are  constants  determined  by  the  method  of  least 
squares  from  the  calibration  data.  In  our  case,  we  substituted  for  R the 
resistance  of  the  thermistor  together  with  the  pad  resistors.  The  errors 
introduced  because  of  this  substitution  are  expected  to  be  small  because 
the  calibration  and  the  use  of  the  chain  was  limited  to  a range  of  5 to 
15°C. 

The  pressure  transducers  together  with  their  bridge/amplifiers  were 
calibrated  in  a pressure  bomb.  There  was  evidence  of  changes  in  the 
calibration  throughout  the  experiment  of  as  much  as  one  meter  in  depth. 
The  chain  was  very  nearly  vertical  while  under  tow  so  that  we  were  able 
to  independently  determine  the  mean  depth  of  the  sensors  during  the 
experiment  from  the  length  of  the  chain  deployed.  The  pressure  sensors 
were  relied  on  to  determine  variations  in  depth. 
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3.  Observations 

Observations  were  made  with  the  towed  thermistor  chain  during  the 
Mixed  Layer  Experiment  (MILE).  The  objective  of  the  MILE  was  to  inves- 
tigate the  response  of  the  upper  ocean  to  atmospheric  forcing.  The 
experiment  was  conducted  in  the  vicinity  of  Ocean  Station  P (SON,  145W) 
during  August  and  September,  1977.  The  experimental  array  is  shown  in 
Fig.  3.  Observations  of  temperature  and  velocity  in  the  upper  100  m were 
made  at  moorings  Ml  and  M2.  Wave  observations  by  use  of  a moored  wave- 
rider  buoy  were  made  at  M3.  In  addition  to  measurements  with  the  towed 
thermistor  chain,  the  shipborne  observations  included:  1)  vertical  pro- 
files of  temperature  and  conductivity  microstructure  (Dillon  and  Caldwell, 
1978;  Dillon  and  Caldwell,  1979;  Caldwell  et  al. , 1979;  2)  towed  measure- 
ments of  temperature,  velocity  and  conductivity  microstructure;  3)  ver- 
tical profiles  of  velocity  microstructure;  4)  towed  saw-tooth  profiles 
of  temperature  and  conductivity;  and  5)  conventional  CTD  profiles.  Ver- 
tical temperature  profiles  were  also  measured  by  use  of  expendable  bathy- 
thermographs dropped  from  aircraft. 

The  thermistor  chain  was  usually  towed  once  around  the  20-km  square 
(Fig.  3)  with  the  midpoint  of  the  instrumented  section  of  the  chain  near 
the  base  of  the  mixed  layer.  The  times  and  locations  of  the  tows  are 
given  in  Table  2 and  Fig.  4.  The  tow  on  5 September  proceeded  four  times 
around  the  five-km  square  centered  on  the  mooring  Ml.  The  tow  tracks 
and  speeds  were  determined  from  satellite  fixes,  radar  fixes  on  the 
moorings  and  dead  reckoning  by  use  of  the  output  from  the  ship's  log  and 
gyrocompass  which  were  recorded  on  our  data  acquisition  system.  During 
tows,  the  ship's  heading  and  speed  were  maintained  as  constant  as 
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possible  on  each  leg.  The  tow  speed  ranged  from  1 to  3 m/s  and  usually 
war  greater  than  2 m/s. 

The  periods  of  the  tows  together  with  wave  height,  surface  tempera- 
ture and  wind  speed  are  shown  in  Fig.  5.  Wind  speeds  during  the  tows 
ranged  from  2 to  12  m/s.  The  highest  wind  speed  during  a tow  occurred 
on  1 September. 

An  example  of  the  temperature  measured  by  use  of  the  thermistor 
chain  is  shown  in  Fig.  6.  The  temperature  oscillations  having  periods 
of  about  10  seconds  are  caused  by  the  thermistor  chain  moving  up  and  down 
in  response  to  the  pitch  and  heave  of  the  ship.  For  a given  amplitude 
of  vertical  motion,  the  amplitude  of  the  temperature  oscillations  is 
proportional  to  the  magnitude  of  the  vertical  temperature  gradient  near 
each  thermistor.  The  upper  thermistors  are  occasionally  in  the  mixed 
layer  where  no  oscillations  are  apparent.  Occasional  patches  of  high- 
frequency  fluctuations,  perhaps  turbulence,  can  be  seen,  e.g.  in  the 
record  from  the  lowest  thermistor  between  2313  and  2314  GMT.  Also  evi- 
dent in  the  records  are  larger-scale  variations  in  temperature  which  are 
taken  to  be  signatures  of  internal  waves.  The  change  in  temperature 
along  the  horizontal  is  occasionally  very  abrupt  as  is  seen  in  the  record 
from  the  lowest  thermistor  between  2304  and  2305  GMT  where  there  is  a 
change  of  nearly  2°C  over  a horizontal  distance  of  30  m. 

The  average  buoyancy  frequency  profile  during  the  experiment  is 
shown  in  Fig.  7.  There  is  a maximum  of  18  cph  in  N at  a depth  between 
25  and  50  m.  Below  this  maximum  there  is  a minimum  in  N of  about  3 cph 
at  a depth  of  approximately  75  m.  A secondary  maximum  in  N of  about  7 
cph  occurs  at  a depth  of  125  m.  Below  200  m,  N is  nearly  constant, 
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about  2 cph.  The  maximum  in  N below  the  base  of  the  mixed  layer  was  much 
larger  than  18  cph  during  and  following  periods  of  high  winds  when 
deepening  of  the  well-mixed  layer  occurred. 

Average  buoyancy  frequency  profiles  over  the  depth  range  of  the 
temperature  measurements  for  each  of  the  tows  are  shown  in  Fig.  8.  Den- 
sity was  determined  from  temperature  by  use  of  a linear  regression  be- 
tween temperature  and  density  observations  between  20  and  50  m depth. 

The  observations  were  taken  by  use  of  a conventional  CTD  operated  from 
the  NOAA  ship  OCEANOGRAPHER.  Observations  were  chosen  which  preceded, 
followed  or  occurred  during  tows  of  the  thermistor  chain.  The  standard 
deviation  of  the  observed  density  about  the  regression  line  was  0.06  ot 
units.  The  errors  in  N introduced  by  not  having  direct  observations  of 
salinity  are  not  expected  to  exceed  10%.  The  largest  values  of  N oc- 
curred on  24  August  following  the  most  intense  storm.  Instantaneous 
profiles  of  N would  show  larger  values  than  those  shown  in  Fig.  8 because 
of  horizontal  variability  along  the  tow  track,  primarily  vertical  dis- 
placements due  to  internal  waves. 

In  most  respects,  the  thermistor  chain  system  performed  well  during 
the  experiment.  Launching  and  retrieving  were  easily  accomplished,  and 
the  fairing  reduced  drag  effectively.  The  winch  and  sheave  performed 
well,  and  bridges,  amplifiers,  and  the  digital  recorder  were  reliable, 
had  low  noise  and  stable  calibrations.  The  biggest  disappointment  was 
the  failure  of  many  of  the  temperature  sensors.  The  failures  were  caused 
by  saltwater  leaks  into  the  potting  between  the  thermistors  and  the  pad 
resistors  and  by  leaks  through  the  glass  coatings  of  the  thermistors. 
These  failures,  which  have  since  been  corrected,  resulted  in  having  only 
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4 to  10  thermistors  working  simultaneously  during  any  tow.  The  vertical 
depth  range  and  density  of  the  measurements  were  therefore  reduced. 
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4.  Analysis 

The  first  step  in  the  analysis  was  to  low-pass  filter  the  time 
series  to  eliminate  the  temperature  oscillations  caused  by  pitch  and 
heave  of  the  ship.  This  was  done  by  computing  sequential,  nonoverlapping 
averages,  each  30  s long.  An  example  of  the  filtered  series  is  shown  in 
Fig.  9.  Isotherm  depths,  usually  at  0.5°C  intervals,  were  then  computed 
by  linear  interpolation  between  the  temperature  records.  Separate  series 
of  isotherm  depths  were  generated  for  each  tow  leg  of  constant  direction. 
Records  of  isotherm  depth  often  could  not  be  computed  for  an  entire  leg 
because  of  the  limited  depth  range  of  the  thermistors  and  the  large  am- 
plitude of  the  internal  waves.  Records  which  were  more  than  80%  complete 
were  completed  by  extrapolation  from  adjacent  isotherms.  This  procedure 
is  justified  by  the  high  coherence  between  adjacent  isotherms  for  scales 
of  1 km  and  larger.  The  resulting  records  varied  in  length  from  50  to 
500  points,  each  point  representing  an  average  over  30  s.  Isotherms  com- 
puted from  the  data  in  Fig.  9 are  shown  in  Fig.  10. 

The  records  of  Isotherm  depth  were  spectrally  analyzed  by  the  use  of 
standard  techniques.  The  Fourier  coefficients  were  computed  directly 
rather  than  by  use  of  a Fast  Fourier  Transform  (FFT)  algorithm.  This  was 
done  because  the  number  of  points  in  a record  was  seldom  near  a power  of 
two  as  required  for  use  of  the  FFT  and  to  avoid  discarding  data  or  intro- 
ducing errors  by  extending  the  records  with  zeros.  Various  windowing, 
detrending  and  prewhitening  methods  were  tried  to  determine  their  effects 
on  the  spectal  estimates.  Frankignoul  (1974)  found  that  the  end  point 
jump  was  an  important  source  of  error  in  spectra  of  internal  waves  com- 
puted from  short  records.  On  the  basis  of  Frankignoul 's  recommendation 


and  our  own  tests,  we  prewhitenend  each  series  by  taking  the  first  for- 
ward difference,  computed  spectra,  and  then  recolored  the  spectrum  by 
dividing  by  the  transform  of  the  differencing  scheme  (Bath,  1974).  We 
did  not  detrend  or  apply  a window  to  the  series.  The  spectra  were 
smoothed  by  averaging  in  non-overlapping  wavenumber  bands,  equally 
spaced  on  a logarithmic  scale.  Frequency  spectra  were  converted  to 
wavenumber  spectra  by  the  use  of  Taylor's  hypothesis  taking  the  mean 
tow  speed  during  a leg  as  the  relevant  velocity. 
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5.  Results 

Spectra  from  20  August  of  the  depth  of  the  10.0  and  11.5°C  Iso- 
therms together  with  the  spectrum  of  depth  from  a pressure  sensor  on 
the  chain  are  shown  In  Fig.  11.  The  spectral  level  of  the  depth  of  the 
chain  is  well  below  the  other  spectra  showing  that  spectra  of  Isotherm 
depth  are  not  signf iciantly  contaminated  by  vertical  motion  of  the 
chain.  The  spectra  of  Isotherm  depth  show  features  typical  of  spectra 
from  other  tows.  The  energy  density  decreases  with  increasing  wave- 
number  except  for  a shoulder  at  a wavelength  of  about  1 km.  The  energy 

_3 

density  is  approximately  proportional  to  k , where  k is  wavenumber,  for 
wavenumbers  greater  than  those  at  the  shoulder.  The  shoulder  does  not 
occur  in  the  spectra  from  all  of  the  tows,  but  there  is  nearly  always 
a change  in  slope  at  a wavelength  of  about  1 km.  The  slopes  of  the 
spectra  in  Fig.  11  together  with  those  from  other  sources  are  summarized 
in  Table  3.  The  spectral  energy  of  the  depth  of  the  10.0°C  isotherm  is 
nearly  an  order  of  magnitude  greater  than  the  spectrum  of  the  U.5°C 
Isotherm  at  all  wavenumbers.  The  11.5°C  isotherm  was  in  a region  of 
weak  vertical  stratification  lying  above  a layer  of  stronger  stratifi- 
cation as  can  be  seen  in  Fig.  10. 

Spectra  of  the  depth  of  the  lowest  and  highest  isotherm  from 
individiaul  tows  were  averaged  to  obtain  the  spectra  shown  in  Fig.  12. 

The  spectra  show  features  similar  to  those  in  Fig.  11,  al through  the 
shoulder  occurring  at  a wavelength  of  1 km  is  not  as  pronounced.  This 
might  be  epxected  since  the  effects  of  doppler-shifting  and  sampling 
packets  of  waves  at  various  angles  to  their  crests  would  tend  to  smear 
the  averaged  spectrum.  The  spectral  energy  of  the  vertical  displacements 


of  the  highest  isotherms  is  lower  than  the  corresponding  energy  of  the 
lowest  Isotherms.  The  difference  is  in  the  same  sense  as  for  the  tow 
on  20  August  (Fig.  11),  although  not  nearly  as  large  and  perhaps  not 
significant. 

The  averaged  spectra  in  Fig.  13  are  scaled  by  the  buoyancy  frequency 
ratio,  N/Nq,  where  N Is  the  local  value  and  NQ  is  taken  to  be  3 cph. 

This  scaling  does  not  collapse  the  averaged  spectra  of  the  depth  of  the 
lowest  and  highest  Isotherms  because,  as  discussed  above,  the  lowest 
isotherm  displacements  tend  to  be  more  energetic  and  in  layers  having 
higher  values  of  N than  the  uppermost  Isotherms  (Fig.  8).  Nevertheless, 
the  two  averaged  spectra  in  Fig.  13  are  in  good  agreement  with  the 
Garrett-Munk  model  (1972,  1975)  despite  the  rapid  vertical  variation  of 
the  local  value  of  N.  There  is,  of  course,  no  break  in  slope  in  the 
model  as  is  suggested  by  the  observations. 

The  averaged  spectra  are  compared  with  towed  spectra  from  other 
sources  in  Fig.  14.  There  is  fair  agreement  between  the  different  sets 
of  observations,  with  spectral  levels  within  a factor  of  five.  The 
spectral  shapes  are  also  similar  with  evidence  of  a break  in  slope  at 
a wavelength  of  about  1 km.  The  slopes  of  the  various  spectra  and 
slope-breaks  are  tabulated  in  Table  3.  Spectra  reported  by  Nelson  and 
Milder  (1979)  from  the  upper  thermocline  also  show  evidence  of  a 
shoulder  at  a wavelength  of  about  1 km.  Pinkel  (1975)  reported  obser- 
vations from  FLIP  in  the  upper  200  m in  which  he  too  found  evidence  of 
a should  In  the  horizontal  wavenumber  spectrum. 

The  spectra  were  examined  for  evidence  of  anisotropy  by  averaging 
spectra  from  different  tow  directions.  Only  the  lowest  isotherm  from 
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each  leg  was  used.  The  results,  plotted  in  Fig.  15,  show  no  evidence 
of  significant  anisotropy.  This  is  in  agreement  with  most  other  obser- 
vations in  the  open  ocean  and  with  the  assumption  of  Garrett  and  Munk 
(1972,  1975).  Nelson  and  Milder  (1979)  did,  however,  find  evidence  of 
anisotropy  at  wavelengths  greater  than  1 km. 

The  spectra  were  also  averaged  according  to  the  local  value  of  N, 
as  shown  in  Figure  16.  Owing  to  the  correlation  between  N and  depth 
(Fig.  8),  the  results  are  similar  to  difference  between  averaged  spec- 
tra of  the  depth  of  the  lowest  and  highest  isotherms  (Fig.  12).  There 
is  evidence  that  for  wavelengths  greater  than  1 km,  higher  values  of  N 
are  associated  with  higher  spectral  energies.  The  usual  N-scaling 
results  in  a greater  spread  of  the  spectral  estimates. 

It  has  been  suggested  (e.g.  Kase  and  Clarke,  1978;  Wunsch,  1975) 
that  internal  waves  may  be  generated  by  wind.  As  a possible  though 
not  necessary  consequence,  one  might  expect  that  spectral  energy  levels 
would  be  higher  during  periods  of  high  winds.  The  spectra  of  the  depth 
of  the  uppermost  isotherm  from  each  tow  were  therefore  grouped  and 
averaged  according  to  the  wind  speed  range  during  and  preceding  each 
tow.  The  results  are  shown  in  Fig.  17.  Spectra  from  tows  during  the 
lowest  wind  speeds  (3-5  m/s)  are  on  average  lower  than  from  tows  during 
higher  wind  speeds.  The  difference  is  not  significant  at  the  95%  level 
and  may  have  been  caused  by  other  factors.  The  spectral  levels  from 
the  two  highest  wind  speed  ranges  are  remarkably  similar. 

Coherence  spectra  between  isotherms  separated  in  the  vertical  were 
computed  from  pairs  of  isotherms  from  each  leg.  Pairs  of  isotherms  which 
were  computed  from  the  same  pair  of  thermistors  more  than  23%  of  the 
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time  were  excluded  from  the  analysis  so  as  to  greatly  reduce  the  possi- 
bility of  introducing  an  artificially  high  coherence.  The  results  from 
the  tow  on  20  August  are  shown  in  Fig.  18.  All  of  the  spectra  are 
qualitatively  similar,  having  high  values  of  coherence  at  the  longest 
wavelengths,  decreasing  with  increasing  wavenumber  and  then  increasing 
again  to  a peak  in  coherence  centered  around  a wavelength  of  1 km.  The 
coherence  between  the  uppermost  (11.5°C)  isotherm  and  the  remaining 
isotherms  is  the  lowest.  The  uppermost  isotherm  was  in  a layer  of  weak 
stratification  overlying  a layer  of  stronger,  nearly  uniform  stratifi- 
cation containing  the  remaining  isotherms  (Fig.  10).  Phase  spectra  be- 
tween pairs  of  isotherms  are  insignificantly  different  from  zero. 

Coherence  spectra  from  individual  tows  were  grouped  and  averaged 
according  to  the  mean  vertical  separation  between  the  isotherms.  The 
results  plotted  in  Fig.  19  are  qualitatively  similar  in  shape  to  Fig.  18. 
With  the  exception  of  isotherm  pairs  separated  by  1 to  3 m,  all  of  the 
spectra  show  a peak  centered  at  a wavelength  of  1 km.  For  separations 
between  1 and  3 m,  coherence  is  high  throughout  the  low  wavenumber  band 
of  the  spectrum,  dropping  only  beyond  a wavelength  of  about  0.5  km. 

Katz  and  Briscoe  (1979)  also  found  evidence  at  350  m depth  of  a 
peak  in  towed  vertical  coherence  between  0.7  and  2 cpkm.  Towed  vertical 
coherence  from  a depth  of  750  m showed  no  evidence  of  a peak.  The 
Garrett-Munk  model  predicts  coherence  decreasing  with  increasing  wave- 
number  in  agreement  with  Katz  and  Briscoe's  (1979)  observations  at  750 
m depth.  Nelson  and  Milder  (1979)  found  higher  values  of  towed  verti- 
cal coherence  in  the  upper  thermocline  than  predicted  by  the  Garrett- 
Munk  model . 


The  peak  In  coherence  at  a wavelength  of  1 km  suggests  that  the 
Internal  wave  spectrum  at  this  wavelength  is  dominated  by  the  low  modes. 
If  higher  modes  are  present,  they  would  result  in  destructive  interfer- 
ence, thereby  decreasing  the  coherence.  This  result  is  consistent  with 
the  analysis  of  fluctuations  of  the  vertical  velocity  of  internal  waves 
by  Dillon  and  Caldwell  (1979)  in  the  upper  200  in  during  MILE.  They 
found  that  the  fluctuations  in  a frequency-band  centered  at  3 cph  were 
dominated  by  the  first  and  second  modes.  The  frequency  spectrum  of 
isotherm  displacement  during  MILE  has  a shoulder  at  about  3 cph  and  then 
falls  off  rapidly  with  increasing  frequency  (Dillon  and  Caldwell,  1979). 
The  frequency  band  centered  at  3 cph  corresponds  to  a wavenumber  band  at 
about  1 cpkm  as  deduced  from  the  dispersion  relation.  Pinkel  (1975)  also 
concluded  on  the  basis  of  observations  in  the  upper  400  m that  the  lowest 
modes  were  dominant  in  the  frequency  band  from  3 to  6 cph. 

We  therefore  tentatively  interpret  our  observations  of  a spectral 
shoulder  and  towed  vertical  coherence  peak  at  a wavelength  of  1 km  as 
evidence  of  the  penetration  from  below  of  low-mode  internal  waves  having 
a frequency  of  about  3 cph. 


6.  Summary 

On  the  basis  of  the  analysis  of  the  vertical  displacement  of  iso- 
therms interpolated  from  towed  thermistor  chain  measurements  between 
20  and  40  m depth  we  can  draw  the  following  conclusions: 

1)  On  average,  the  spectrum  of  isotherm  displacement  decreases 

-3/2 

with  increasing  wavenumber  proportional  to  k out  to  1 cpkm. 
Above  1 cpkm  the  spectral  slope  is  between  -2  and  -3. 

2)  There  is  often  a shoulder  in  the  spectra  of  isotherm  displace- 
ment at  a wavelength  of  1 km.  This  feature  is  smeared  in  the 
spectra  averaged  from  many  tows. 

3)  There  is  no  evidence  of  significant  horizontal  anisotropy  in 
spectra  averaged  over  many  tows. 

4)  Averaged  spectra  of  isotherm  displacement  appear  to  be  inde- 
pendent of  the  local  buoyancy  frequency.  Scaling  individual 
spectra  by  N,  as  is  often  done,  therefore  increases  the 
scatter. 

5)  The  spectra  of  isotherm  displacement  does  not  depend  signi- 
ficantly on  wind  speed  for  wind  speeds  below  10  m/s. 

6)  Despite  the  rapid  vertical  variation  of  N,  the  spectra  of 
isotherm  displacements  are  in  fair  agreement  with  the 
Garrett-Munk  model  (1972,  1975)  and  with  spectra  reported 
by  other  investigators. 

7)  There  is  a peak  at  a wavelength  of  1 km  in  coherence 
spectra  between  pairs  of  isotherms  separated  in  the  vertical. 
The  peak  is  interpreted  as  evidence  that  the  wavenumber  band 
at  1 cpkm  is  dominated  by  the  low  modes. 
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TABLE  1 

. Specifications  of  Towed  Chain  Measurements 
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Sensor 

Manufacturer  Accuracy 

Resolution 

Response 

Time 

Thermistor, 
glass-coated  bead, 
0.2  mm  diameter 

Thermometries  4 0.01°C 

Model  P 85  BIO 

± o.oorc 

0.12s 

Pressure 

(Depth) 

Kulite 

Model  PTQH-360-250  ± 0.5  m 

4 0.05  m 

Is 
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TABLE  2.  Summary  of  Towed  Thermistor  Chain  Observations 


Date 

Time 

Begin 

(GMT) 

Duration 

(Min) 

Tow  Speed 
(m/s) 

Track 
(see  Fig 

20  Aug  77 

1846 

195 

1.6 

A 

24  Aug  77 

1654 

465 

1. 6-2.0 

B 

26  Aug  77 

1521 

534 

2. 3-3.0 

C 

28  Aug  77 

1515 

718 

1.1-3. 3 

C 

30  Aug  77 

1736 

575 

2. 5-3. 2 

C 

1 Sep  77 

1723 

440 

2. 2-2. 3 

D 

3 Sep  77 

1811 

520 

2. 3-3.1 

C 

5 Sep  77 

1817 

487 

2.6 

E 

8 Sep  77 

1814 

173 

2.8 

F 

TABLE  3.  Spectral  Slopes 


SOURCE 

SLOPE 

SLOPE  BREAK 

Low  k 

High  k 

(c/km) 

Run  1,  10° 

-1.3 

-3.0 

1 .0 

Run  1,  11.5° 

-1.5 

-3.6 

1.0 

MILE,  Ensemble 

-1.6 

-2.2 

0.8 

Katz  (1973) 

-1.5 

-2.3 

1.8 

Nelson  and  Milder 

-1.5 

-2.5 

1.3 

(1979) 

sale SPafciJm 
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FIGURE  LEGENDS 

1.  Schematic  diagram  of  the  thermistor  chain  under  tow. 

2.  A diagram  showing  the  installation  of  a thermistor  and  bridge/ampli- 
fier in  a modified  section  of  fairing.  The  dimensions  of  the  fairing 
are  15.5  x 2 x 10  cm.  The  nose-pieces  fit  over  a steel  cable  and 
overlap  with  adjacent  tail  sections. 

3.  The  MILE  array.  Station  P (50N,  145W)  is  at  northern  apex.  Moorings 
are  at  Ml,  M2,  and  M3.  The  thermistor  chain  was  towed  by  the  NOAA 
ship  OCEANOGRAPHER,  usually  around  the  20  km  square. 

4.  Idealized  thermistor  chain  tow  tracks.  Clockwise  tows  were  made 
about  the  squares,  generally  beginning  near  the  southwest  corner. 

Tow  A was  made  to  the  northwest,  F to  the  southwest.  See  Table  2 
for  time  and  towing  speed. 

5.  Near-surface  observations  during  the  MILE.  The  temperatures  at  5 m 
depth  are  hourly  averages  from  the  Ml  mooring.  The  wind  speeds  are 
hourly  averages  from  the  OCEANOGRAPHER.  Wave  heights  are  from  the 
wave-rider  at  M3, 

6.  Time  series  of  thermistor  outputs  recorded  on  3 September  1977.  The 
data  were  recorded  at  8 samples/s  and  plotted  over  0.5  s.  The  ship 
was  proceeding  southward  at  3.0  m/s. 

7.  Buoyancy  frequency  vs.  depth  averaged  over  the  entire  experiment. 

8.  Buoyancy  frequency  profiles  determined  from  the  thermistor  chain 
data  averaged  over  entire  tows.  Numbers  on  the  plot  refer  to  the 
date  of  tow,  two  digits  referring  to  August  1977  and  one  digit  to 
September. 
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9.  Low  pass  filtered  temperature  time  series  of  the  tow  on  20  August 
1977.  The  ship  was  traveling  to  the  northeast  at  1.6  m/s. 

10.  Isotherm  depths,  20  August  1977,  computed  by  linear  Interpolation 
to  the  data  in  Fig.  9. 

11.  Spectra  of  the  depth  of  the  10.0°  and  11.5°  isotherms  and  the  spectrum 
of  depth  from  a pressure  sensor  on  the  chain.  The  tow  occurred  on 

20  August.  Average  depths  were  37.2,  27.3,  and  33.8  m,  for  the 
isotherms  and  pressure  sensor  respectively.  Symbols  are:  +,  10.0°C 
isotherm;  X,  11.5°C  isotherm;  and  □,  pressure.  Vertical  bars  repre- 
sent the  95%  confidence  intervals. 

12.  Ensemble- averaged  spectra  of  the  depth  of  the  lowest  (+)  and  highest 
isotherms  (X)  for  the  entire  experiment.  Averaging  was  performed 
without  regard  to  actual  isotherm  values  or  depths.  Vertical  bars 
represent  the  95%  confidence  intervals. 

13.  Comparison  of  ensemble-averaged  N-scaled  spectra  of  the  depth  of  the 
lowest  and  highest  isotherms  and  the  Garrett-Munk  (1972,  1975)  models. 

y 

Local  values  of  N were  derived  from  averaged  thermistor  data  for 
each  time  series.  Nq  is  3 cph. 

14.  Comparison  of  N-scaled  spectra  of  isotherm  depth  from  several  sources. 

15.  Ensemble-averaged  spectra  of  isotherm  depth  from  tows  in  various 
directions.  Only  the  lowest  isotherm  from  each  tow  was  used.  Sym- 
bols are:  X,  north;  +,  east;  □,  south;  and  O,  west.  Vertical 
bars  represent  the  95%  confidence  intervals. 

16.  Ensemble-averaged  spectra  of  isotherm  depth  grouped  according  to  the 
local  value  of  N.  Symbols  are:  X,  N < 5;  +,  5 < N < 10;  a, 

10  < N < 15;  O.  N > 15  (all  units  in  cph).  Vertical  bars  represent 
95%  confidence  limits. 
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17.  Ensemble-averaged  spectra  of  the  depth  of  the  highest  Isotherms  for 
tows  grouped  according  to  average  wind  speed.  Symbols  are;  X, 

8-10  m/s;  +,  5-7  m/s;  and  □,  3-5  m/s  Vertical  bars  represent 

95%  confidence  limits. 

18.  Coherence  spectra  between  Isotherms  from  the  tow  on  20  August.  The 
solid  line  Is  the  95%  significance  level. 

19.  Ensemble-averaged  coherence  spectra  between  isotherms  having  various 
mean  vertical  separations,  regardless  of  absolute  depth.  The  solid 
line  is  the  95%  signficance  level. 
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APPENDIX  A 
MILE  Tow  Parameters 

This  section  contains  listings  of  parameters  associated  with  the 
thermistor  chain  tows  and  some  elementary  results.  During  the  experiment, 
tows  were  made  on  a total  of  10  occasions.  Data  from  the  sensors  were 
recorded  In  continuous  sections  called  runs,  with  a single  tow  comprising 
of  one  or  more  runs.  The  10  tows  were  recorded  in  a total  of  24  runs. 

For  analysis,  each  tow  was  broken  up  into  legs,  which  are  tow  segments 
of  nearly  constant  direction  and  speed.  Run  and  leg  boundaries  had  no 
specific  relationship  so  that  a single  leg  could  span  one  or  more  runs, 
or  a single  run  could  comprise  one  or  more  tow  legs.  The  numbering  of 
the  legs  Is  derived  from  the  run  numbers,  and  is  best  described  by  exam- 
ple. Run  4 and  part  of  run  5 make  up  leg  4-5,  leg  15A  is  contained  in 
the  first  half  of  run  15  and  so  on.  A total  of  37  tow  legs  were  recorded. 

One  of  the  parameters  obtainable  from  the  thermistor  chain  data  is 
the  buoyancy  frequency,  N.  This  is  the  natural  frequency  of  the  vertical 
oscillation  of  a water  parcel  in  a stably  stratified  environment  and  is 
given  by 

N2  = . a Ifi. 

p az 

where  g is  the  acceleration  due  to  gravity,  p is  the  density  and  z is 
the  vertical  space  coordinate,  positive  upwards.  Since  conductivity  was 
not  measured  on  the  chain,  a relation  between  temperature  and  density 
was  needed.  This  was  obtained  from  CTD  profiles  measured  on  the  NOAA 
ship  OCEANOGRAPHER  by  the  group  led  by  Stan  Hayes  of  the  Pacific  Marine 
Environmental  Laboratory.  Profiles  were  selected  which  were  taken  during 


or  within  a few  hours  oi  the  thermistor  chain  tows.  Values  of  tempera- 
tu.e,  salinity,  and  from  every  10  m in  the  depth  range  -20  to  -50 
were  assembled.  A linear  relation  between  temperature  and  was  derived 

ot  = -0.19179  T (°C)  + 27.038 

The  rms  deviation  of  the  observations  from  this  line  was  0.0647  units. 
Using  this  relation  and  the  horizontally  averaged  temperature  of  each 
thermistor  for  each  tow  leg,  profiles  of  were  calculated.  Using  the 
relation  p = 1.0  + 1000.  x o^,  the  values  of  N at  the  midpoints  between 
thermistors  were  obtainable.  Just  as  the  isotherm  depths  were  determined 
by  linear  interpolation  of  the  thermistor  data,  the  average  values  of  N 
for  each  isotherm  were  found  by  linear  interpolation  of  the  values  at  the 
midpoints  between  thermistors. 


MILE  Thermistor  Chain  Tow  Parameters 


Run 

Data 

Tima 

(CMT) 

Tow 

Spaad 

(M/S) 

Tow 

Dlractlon 

(Dag) 

Duration 

Tima  Dlatanca 
(Min)  (ICM) 

Wind 

Spaad 

(M/S) 

Wind 

Dlractlon 

(Dag) 

Wind  Wava 
Haight 
(M) 

1 

20-AUC-77 

1846 

1.6 

NE 

195 

18.5 

5 

120 

0.5 

2 

22-AUC-77 

1648 

1.1 

W 

6 

0.4 

18 

240 

3.0 

3 

24-AUC-77 

1654 

1.6 

SNW 

114 

10.9 

5 

310 

1.5 

4-5 

It 

1852 

1.6 

NW 

171 

16.4 

4 

310 

1.5 

5-6 

It 

2142 

2.0 

S 

180 

21.6 

5 

300 

1.5 

7 

26-AUG-77 

1521 

2.3 

N 

162 

22.  7 

6 

250 

0.5 

8 

it 

1810 

2.7 

E 

103 

16.5 

6 

240 

0.8 

9A 

it 

1959 

3.0 

S 

149 

26.4 

a 

260 

1.0 

98 

•t 

2229 

3.0 

W 

120 

21.2 

7 

260 

1.0 

10-11 

28-AUG-77 

1515 

1.1 

u 

26 

1.7 

8 

250 

1.0 

11A 

" 

1541 

2.7 

N 

123 

20.2 

7 

250 

1.5 

UB 

it 

1744 

3.3 

E 

81 

15.9 

7 

280 

1.5 

11-12 

" 

1905 

1.9 

S 

240 

27.7 

8 

250 

1.5 

12A 

•t 

2111 

2.5 

S 

112 

16.7 

8 

240 

1.5 

128 

it 

2303 

2.0 

w 

136 

16.3 

7 

240 

1.5 

13 

30-AUG-77 

1736 

3.0 

s 

27 

4.9 

8 

060 

1.0 

13-14 

it 

1803 

2.5 

N 

158 

23.6 

8 

050 

1.0 

14 

it 

2041 

3.2 

E 

120 

22.7 

8 

040 

1.0 

15A 

•t 

2258 

3.2 

S 

150 

28.8 

7 

040 

1.0 

158 

31- AUG- 7 7 

0128 

3.0 

u 

120 

21.6 

5 

040 

1.5 

16 

I-SEP-77 

1723 

2.3 

N 

25 

3.4 

12 

360 

2.5 

17 

•t 

1757 

2.2 

N 

190 

25.2 

10 

350 

2.5 

17-18 

2107 

2.2 

E 

225 

29.7 

8 

340 

2.0 

19 

3-SEP-77 

1811 

2.3 

N 

130 

17.9 

7 

360 

0.5 

20 

•• 

2025 

2.6 

E 

118 

18.6 

4 

340 

0.5 

21A 

tt 

2227 

3.1 

S 

140 

25.6 

5 

330 

0.5 

218 

4-SEP-77 

0047 

2.9 

W 

132 

22.6 

6 

320 

0.5 

22A 

5-SEP-77 

1817 

2.6 

u 

16 

2.5 

3 

280 

0.5 

22B 

•t 

1833 

2.6 

N 

33 

5.2 

3 

270 

0.5 

22C 

it 

1906 

2.6 

E 

21 

3.3 

3 

270 

0.5 

2 2D 

•t 

1927 

2.6 

S 

27 

4.2 

3 

280 

0.5 

22E 

it 

1954 

2.6 

V 

33 

5.2 

3 

290 

0.5 

22F 

•t 

2027 

2.6 

N 

37 

5.8 

3 

300 

0.5 

22C 

it 

2104 

2.6 

E 

22 

3.4 

3 

300 

0.5 

22H 

it 

2126 

2.6 

S 

29 

4.5 

3 

300 

0.5 

221 

tt 

2155 

2.6 

W 

30 

4.7 

3 

300 

0.5 

22J 

it 

2225 

2.6 

N 

32 

5.0 

3 

300 

0.5 

22K 

it 

2257 

2.6 

E 

28 

4.4 

2 

300 

0.5 

22L 

tt 

2325 

2.6 

S 

29 

4.5 

3 

290 

0.5 

22M 

it 

2354 

2.6 

w 

32 

5.0 

3 

290 

0.5 

23A 

6-SEP-77 

0031 

2.6 

N 

28 

4.4 

3 

300 

0.5 

238 

it 

0059 

2.6 

E 

29 

4.5 

3 

300 

0.5 

23C 

it 

0128 

2.6 

S 

31 

4.8 

3 

310 

0.5 

230 

tt 

0159 

2.6 

w 

30 

4.7 

2 

310 

0.5 

24 

8-SEP-77 

1814 

2.8 

sw 

173 

28.6 

5 

250 

1.5 
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MILE  Depths  (a)  of  Working  Theralstors 


CHANNEL 


Run  IS 

16 

18 

19 

20 

21 

22 

23 

24 

27 

30 

1 

22.74 

23.79 

24.84 

25.89 

26.93 

27.98 

29.03 

32.17 

35.32* 

2 

21.69* 

22.74 

24.84 

25.89 

26.93 

27.98 

29.03 

30.08 

3 

24.84 

25.89 

26.93 

27.98 

30.08 

4 

24.84 

25.89 

26.93 

27.98 

29.03 

30.08 

S 

24.84 

25.89 

26.93 

27.98 

29.03* 

30.08 

6 

24.84 

25.89 

26.93 

27.98 

30.08 

7 

21.69 

22.74 

24.84 

25.89 

26.93 

27.98 

30.08 

8 

21.69* 

22.74 

24.34 

25.89 

26.93 

27.98 

30.08 

9 

22.74 

24.84 

25.89 

26.93 

27.98 

30.08 

10 

24.84 

25.89 

26.93 

27.98 

30.08 

11 

22.74 

24.84 

25.89 

26.93 

27.98 

30.08 

12 

22.74 

24.84 

25.89 

26.93 

27.98 

30.08* 

34.27 

13 

28.51 

30.60 

31.65 

32.70 

33.75 

40.03 

14 

28.51 

30.60 

31.65 

32.70 

33.75 

40.03 

15 

28.51 

30.60 

31.65 

32.70 

33.75 

40.03 

16 

18.55 

20.65 

22.74 

23.79 

30.08 

17 

31.86 

33.96 

36.05 

37.10 

43.39 

18 

31.86 

33.96 

36.05 

37.10 

43.39 

19 

29.13 

31.23 

33.33 

34.37 

20 

29.13 

31.23 

33.33 

34.37 

37.52 

21 

29.13 

31.23 

33.33 

34.37 

37.52* 

22 

28.56 

30.65 

32.75 

33.80 

40.09 

23 

27.51* 

28.56 

30.65 

32.75 

33.80 

40.09 

24 

29.66 

31.75 

33.85 

34.90 

41.19 

* 

Working  for  part 

of  run  only 

Run 

Chan 

Working 

30  Sac  Points 

Total 

30  Sac  Points 

1 

30 

1 

259 

391 

2 

15 

1 

7 

13 

5 

22 

1 

111 

310 

8 

15 

1 

98 

206 

12 

23 

1 

440 

497 

21 

24 

1 

457 

545 

23 

15 

55 

237 

237 

50 


RUN 

1 


■> 


1 1 A 


11B 


12A 


AVERAGE 
Z (M) 


-23.27 
-24.31 
-25.36 
-26.41 
-27.46 
-28.51 
-30 . 60 
-33.75 
-38.46 


-22.22 
-23.79 
-25.36 
-26.41 
-27 . 46 
-28.51 
-29 . 55 


-25.36 

-26.41 

-27.46 

-28.51 

-29.55 


-23 . 79 
-25 . 36 

-26.41 

-27.46 

-29.03 


-23.79 
-25.36 
-26.41 
-27.46 
-29 . 03 


-23,79 

-25.36 

-26.41 

-27,46 

-29.03 

-32.17 


AVERAGE 
N (CPH) 


5.08 
5.02 
4 . 89 
5.08 
7.32 
9.72 
8.49 
12.05 
1 2 . 22 


16.06 
10.04 
1 1 . 36 
10.42 
8.87 
8.39 
7.36 


1 .49 
3 . 22 
8.01 
17.47 
20.49 


6.39 
7.42 
9 . 58 
12.15 
14.88 


2.44 

2.70 

3.47 

4.95 

7.99 


5 . 70 
4.59 
3.56 
3.42 
6.45 
18.34 


STD  DEV 
N ( CF'H ) 


1.61 
1 .61 
1.97 
2 . 1 6 
2.03 
2.63 
3.95 
4,12 
2.32 


1 .96 
1.51 
1 .03 
1 .37 
1 . 70 
1.64 
1.64 


2.43 
5 . 0 1 
7.59 
7.75 
4.95 


3.97 

6.31 

6.93 

5 . 72 


1.04 

1,70 

1.72 

2.44 

4.69 


1.51 

2.25 

2.11 

1.66 

3,95 

4.09 


# 30  SEC 
POINTS 


i 

391 

391 

391 

391 

39.1 

391 

259 

259 


7 

13 

13 

13 

13 

13 

13 


270 

270 

180 

180 


246 

246 

246 

246 

246 


163 

163 

163 

163 

163 


224 

224 

224 

224 

224 


51 


RUN 

AVERAGE 

AVERAGE 

STD  DEV 

# 30  SEC 

Z <M) 

N (CPU) 

N (CPH) 

POINTS 

12B 

-23.79 

5.04 

1 .90 

273 

-25 . 36 

4.07 

2.07 

273 

- 26.41 

3 . 68 

2.19 

773 

-27.46 

3 . 88 

2.77 

273 

- 29.03 

7.98 

3.30 

216 

-32.17 

10.82 

4.34 

216 

13 

-29.55 

6.58 

2.89 

54 

31 . 13 

10.82 

8 . 70 

54 

-32.17 

14 . 22 

7.06 

54 

-33.22 

16.54 

5.30 

54 

-36.89 

17.02 

2.21 

54 

13-14 

-29 . 55 

8 . 4 1 

4 . 54 

3 1 0 

-31 . 13 

1 0 . 28 

6.76 

310 

-32.  17 

12.28 

6.21 

310 

-33.22 

13.37 

5.37 

310 

-36.89 

15.92 

3.22 

3 1 0 

14 

-29 . 55 

8 . 39 

4.04 

241 

-31 .13 

1 1 . 38 

5.81 

24  1 

-32.17 

14.24 

6.24 

241 

-33.22 

15.73 

6.13 

241 

- 36 . 89 

15.71 

2.13 

241 

15A 

-29 . 55 

4.14 

2.36 

300 

-31 . 13 

3.93 

3.47 

300 

-32.17 

4.59 

4.04 

300 

-33.22 

5 . 50 

4.87 

300 

-36.89 

14.93 

4.24 

300 

15B 

-29.55 

6.73 

4.34 

241 

-31 . 13 

8.91 

6 . 66 

241 

-32.17 

1 1 . 36 

7.38 

241 

-33 . 22 

13.36 

7.08 

241 

-36.89 

15.98 

2.74 

241 

16 

- 19.60 

l . 33 

0.14 

50 

-21 .69 

0.57 

0,30 

50 

-23.27 

1 . 33 

0.28 

50 

-26.93 

3 . 88 

t .12 

50 

52 


RUN 

AVERAGE 

AVERAGE 

STD  DEV 

# 30  SE 

Z (M) 

N (CPU) 

N ( CPH ) 

POINTS 

17 

-32.91 

16.95 

4.71 

380 

-35.00 

16.09 

4.47 

380 

-36.58 

14.95 

4.31 

380 

-40.24 

12.26 

2.70 

380 

17-18 

-32.91 

9.83 

6.57 

440 

-35.00 

13.22 

6.50 

440 

-36.58 

15.47 

6.13 

440 

-40.24 

15.62 

3.34 

440 

19 

-30.18 

7.26 

7 . 20 

261 

-32.28 

8.41 

6.94 

261 

-33.85 

9.84 

5.97 

261 

20 

-30.18 

8,33 

6.07 

237 

-32.28 

12.81 

6.08 

237 

-33 . 85 

16.10 

6.05 

237 

-35.95 

18.74 

3.14 

237 

21 A 

-30.18 

6.64 

6.51 

280 

-32.28 

10.34 

5.98 

280 

-33 . 85 

14.1 1 

5.50 

280 

-35 . 95 

16 . 72 

3.48 

280 

21 B 

-30. 18 

11.62 

6.45 

265 

-32 . 28 

17.63 

4 . 94 

265 

~33 . 85 

20 . 7 1 

4.73 

265 

-35.95 

18.17 

2.57 

177 

22A 

-29.61 

1.58 

0.73 

32 

-31.70 

1 . 98 

0.85 

32 

-33 . 27 

3 . 1 4 

.1  .54 

32 

-36.94 

9 . 99 

2.16 

32 

22B 

-29.61 

2,41 

0.64 

66 

■31 .70 

3.18 

1 .76 

66 

■33.27 

5.24 

2.35 

66 

■36 . 94 

12,06 

2 . 3 1 

66 

r 


i 


RUN 

AVERAGE 

AVERAGE 

STD  DEV 

* 30  SE 

Z < M) 

N ( CPH ) 

N ( CPH ) 

POINTS 

22K 

-29.61 

5.02 

2.13 

56 

-3 1 . 70 

10.20 

4.33 

56 

-33.27 

13.01 

3.31 

56 

-36.94 

15.61 

1.63 

56 

22L 

-29.61 

3.49 

1.06 

58 

-31 .70 

3.69 

2.99 

58 

-33.27 

5.26 

3.49 

58 

-36.94 

14.46 

3.23 

58 

22M 

i“ 4 

>0 

1 

3.15 

0.75 

65 

-31.70 

3.97 

2.75 

65 

-33.27 

6.54 

3.77 

65 

-36.94 

15.66 

2.13 

65 

23A 

-29.61 

9.27 

6.14 

56 

-31.70 

12.29 

3.48 

56 

-33.27 

12.88 

2.36 

56 

-36.94 

16.58 

0.95 

56 

23B 

-28.03 

9.38 

5.68 

58 

-29.61 

13.76 

3.98 

58 

-31.70 

15.86 

2.94 

58 

• 33.27 

14.95 

3.63 

58 

-36.94 

15.39 

1.94 

58 

23C 

-28.03 

5.87 

1 . 26 

62 

29.61 

10.02 

5.22 

62 

-31 < 70 

17.22 

5.80 

62 

-33.27 

17.66 

3.94 

62 

-36.94 

15.19 

1.79 

62 

23D 

-28 . 03 

13.19 

5.28 

61 

-29.61 

1 4 . 75 

5.29 

61 

-3 1 . 70 

18.07 

3.99 

61 

-33.27 

19.79 

2.98 

61 

-36 . 94 

14.10 

3.58 

6.1 

1 


55 


RUN 

AVERAGE 

AVERAGE 

STD  DEV 

# 30  SEC 

Z (M> 

N (CPU) 

N ( CPH ) 

POINTS 

24 

-30.71 

8.69 

4.31 

34  7 

-32.80 

10.42 

5.29 

34  7 

-34.37 

11.72 

5.41 

347 

--38.04 

15.09 

1 .95 

347 

•V** 


56 


RUN 

# 30  SEC 

ISOTHERM 

AVERAGE 

STD  DEV 

N 

AVERAGES 

(DEG  C) 

Z (M) 

Z < M ) 

(CRH) 

1 

390 

11.50 

-27.28 

1.89 

6.93 

11.00 

-32.34 

2.45 

10.46 

10.50 

-34.94 

3.14 

12.09 

10.00 

-37.17 

3.84 

12.17 

9 

12 

12.00 

-22.98 

0.67 

13.14 

11.50 

-25.50 

1.01 

11.23 

11.00 

-28.29 

1 . 15 

8.49 

4-5 

274 

12.00 

-27.67 

1.25 

9.97 

11.90 

-28.32 

1.37 

1 5 . 78 

HA 

246 

1 1 . 90 

-25.94 

1.62 

8.61 

11.50 

-28.08 

1 .84 

13.23 

11B 

1*2 

12.25 

-26.92 

1 . 82 

4 . 1 9 

12A 

224 

12.00 

-26.35 

1.62 

3 . 6 1 

1 1 . 90 

-29,97 

0.99 

10.01 

1 1 . 50 

-30,80 

1.13 

13.13 

11  .00 

-31.73 

1 . 35 

16.67 

10.50 

-32.57 

1.47 

1 9 . 82 

10.00 

-33.35 

1.53 

22,77 

1 2E< 

272 

11.80 

-29.90 

1 . 85 

8.76 

13 

54 

11.50 

-32.29 

1 . 44 

14.48 

11.00 

-33 , 30 

1 . 53 

1 6 . 55 

10.50 

-34. 11 

1 , 79 

16.66 

10.00 

-34,90 

2.10 

16.76 

9.50 

-35,72 

2.34 

16.87 

9.00 

-36,60  . 

2.50 

16.98 

1 3-  .1 4 

316 

1 1 . 50 

-31 . 98 

2 . 77 

11.81 

1 1 . 00 

-32.99 

2.74 

1 3.13 

1 0 . 50 

-33.89 

2.84 

13.84 

10.00 

-34.84 

2.86 

14 , 49 

9.50 

-35.82 

2 . 88 

15.1 8 

9.00 

-36.90 

2 . 80 

1 5 , 9 '2 

RUN 

* 30  SEC 

ISOTHERM 

AVERAGE 

STD  DEV 

N 

AVERAGES 

(DEG  C) 

Z (M) 

Z < M ) 

(CRH) 

14 

240 

11.50 

-32.13 

2.10 

14.12 

11.00 

-33.27 

2. 13 

15. 73 

10.50 

-34.28 

2.30 

15.72 

10.00 

-35.22 

2.48 

15.72 

9.50 

-36.24 

2.58 

15.71 

9.00 

-37.31 

2.59 

15.71 

15A 

300 

1 1 . SO 

-34.95 

1.80 

9.95 

11.00 

-36.47 

2.51 

13.84 

15B 

240 

11.50 

-32.40 

2.18 

11 .79 

1 1.00 
10.50 

-33.60 

-34.60 

2.06 

2.13 

13.63 

14.35 

10.00 

-35 ,57 

2.28 

15.04 

9.50 

-36.53 

2.44 

15.72 

9.00 

-37.51 

2 . 53 

16.43 

16 

50 

12.25 

-25.68 

1.24 

3.00 

17 

380 

9.00 

-34,78 

2.93 

1 6.18 

8.50 

-35.96 

3.00 

15.39 

8 . 00 

-37.38 

3.07 

14.36 

7.50 

-39.18 

2.85 

13.04 

7 . 00 

-41 . 18 

2.53 

1 1.57 

17-18 

450 

11.50 

-35,33 

3 . 08 

13.69 

11.00 

10.50 

-36.28 

-37.08 

3.11 

3 . 27 

15.04 

15.49 

.10.00 

-37.91 

3.40 

1 5 , 52 

9.50 

-38.78 

3.51 

t 5 .56 

9.00 

-39.70 

3.52 

15.59 

19 

190 

12.35 

-31 .46 

1 . 94 

7 . 96 

.'0 

236 

12.00 

-32.08 

2.09 

.1  2 , 39 

1.1 .50 

-33.32 

2 . 1 1 

1 4 . 99 

1 1.00 

-34.33 

2.16 

1 6 . 7 1 

10.50 

-35.22 

2 . 2 1 

i 7.  a? 

# 30  SEC 

ISOTHERM 

AVERAGE 

STD  DEV 

N 

AVERAGES 

< DEG  C) 

2 (M) 

7 <M> 

(CRH) 

280 

12.00 

-32.55 

2 . 36 

10.99 

1 1 .50 

-33.67 

2.51 

13.69 

11.00 

-34.59 

2.64 

15,04 

264 

11.50 

-31.56 

1.68 

15.56 

11.00 

-32.40 

1.65 

17.87 

10.50 

-33.08 

1.66 

19.19 

10.00 

-33.68 

1.67 

20.38 

9.50 

-34,3.1 

1 . 72 

20.16 

32 

12.35 

-32.92 

1 . 66 

2.88 

12.30 

-34.09 

0.67 

4.67 

12.20 

-34.94 

0.91 

6 . 26 

12.00 

-36.30 

1.51 

8.80 

66 

1 2 . 30 

-32.57 

1 .29 

4.32 

12.20 

-33.88 

0 . 78 

6 . 36 

12.00 

-35.04 

0.99 

8 . 52 

11,50 

-37.22 

1.71 

12.58 

42 

12.35 

-30.11 

1 .42 

3 . 25 

12,30 

-32.29 

1 . 52 

5 . 0 1 

12.20 

-33 . 38 

1 . 20 

6.51 

12.00 

-34.40 

0.99 

8.45 

1 1 . SO 

-36.28 

1 . 15 

12.00 

11.00 

-37.99 

1 .44 

1 5 . 23 

54 

12.35 

-31.88 

1 . 16 

2.33 

12.30 

-33 . 58 

0.51 

5 . 1 1 

12.20 

-34.04 

0.20 

6.62 

12.00 

-34.53 

0.24 

8.24 

11.50 

-35.71 

0.49 

12.16 

11.00 

-36.90 

0.75 

16.08 

10.50 

-38.08 

1 . 00 

19.99 
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RUN 

# 30  SEC 

ISOTHERM 

AVERAGE 

STD  DEV 

N 

AVERAGES 

(DEG  C> 

Z (M) 

Z (M) 

( CPH ) 

22E 

66 

12.35 

-31.23 

1 .59 

4.65 

12.30 

-32.35 

1.36 

6.97 

12.20 

-33.08 

1 . 15 

9.10 

12.00 

-33.81 

0,97 

1C. 62 

11.50 

-35.13 

0.93 

13.03 

11.00 

-36.30 

1 . 16 

15.16 

10.50 

-37.45 

1.36 

1 7 . 26 

10.00 

-38.56 

1.50 

19.27 

22F 

74 

12.30 

-30.66 

1.22 

5,78 

12.20 

-32.10 

1 .28 

7.89 

12.00 

-33.27 

1 . 18 

9.08 

11.50 

-34.91 

1 .01 

12.22 

11 .00 

-36.17 

1.27 

14.62 

10.50 

-37.37 

1 . 55 

16.92 

10.00 

-38.49 

1.71 

19.07 

22G 

44 

12.35 

-31.87 

1.15 

2.34 

12.30 

-33.78 

0.81 

3.86 

12.20 

-34.97 

0.83 

5.77 

12.00 

-36.74 

1.61 

8.60 

22H 

58 

12.35 

-32.90 

1,22 

2.61 

12.30 

-34.16 

0 . 92 

4.87 

12.20 

-35.04 

1 , 64 

6 , 85 

12.00 

-36.09 

2.02 

9 . 1 3 

221 

60 

1 2 . 35 

-34.86 

1.96 

4,14 

22  J 

64 

12.35 

-30.10 

1 . 74 

4.29 

12.30 

-32 . 20 

2.04 

6.18 

1 2 , 20 

-33.22 

1 . 96 

7 . 46 

12.00 

-34.41 

2 . 33 

9,21 

11.50 

-36.54 

3.01 

12.34 

60 


RUN 

• JO  SEC 

ISOTHERM 

AVERAGE 

STD  DEV 

N 

AVERAGES 

< DEG  C> 

/.  <M) 

7 < M ) 

< CPU ) 

22N 

56 

12.  JO 

“JO.  IV 

1.00 

6.47 

12.20 

-31  . *14 

1 . 18 

V , 56 

12.00 

-32. JV 

1 . 22 

11.42 

1 1 . 50 

“33.96 

1 . 24 

1 3 . 50 

1 1 . 00 

“35 . 1 V 

t . 47 

14.37 

1 0 . 50 

36 . 38 

1 , 6 1 

15.21 

1 0 . 00 

-37.62 

1.69 

16.09 

9 . 50 

-38.83 

1,75 

,16.96 

22L 

58 

12.  J 5 

-31 ,38 

2.15 

3 .66 

12.  JO 

-32. /J 

1 . 98 

4.71 

12.20 

-*33.62 

1 .62 

6.13 

12.00 

-34.44 

1 . 55 

8.19 

1 1 , 50 

-36 . 05 

1.65 

12.21 

1 1 . 00 

-37.43 

1.88 

1 5 . 68 

1 0 . 50 

-38 . 7 1 

1.99 

1 8 . 88 

;.*2M 

64 

12.  J 5 

31 .45 

t .54 

3.88 

12,  JO 

32.70 

l . 33 

5 . 60 

12.20 

“33,44 

0 . 98 

6.94 

12.00 

-34.21 

0.  72 

8.85 

1 1 .50 

35 , 58 

0 . 8. 

12.77 

11,00 

■36.88 

1.13 

15.52 

10.50 

38.  t 1 

l .32 

1 8 . 56 

ja 

56 

1 2 . 00 

30.84 

1 .85 

1 1 , 05 

1 1 , 50 

“32*27 

2 . 05 

1 7.50 

t 1 . 00 

-33.52 

7 . 08 

1 3 . 1 3 

1 0 . 50 

■34,72 

2 . 0 7 

1 4 . 34 

1 0 . 00 

-35.86 

2 . 03 

15.49 

V . 50 

37.03 

1 . 92 

t 6 • 6 6 

0.00 

38.09 

1.91 

17.  74 

ojr 

58 

.1  1 , 50 

79 . 98 

1 .58 

14,13 

t 1 . 00 

31 . 75 

1 . 0 

1 5 , 4 t 

10.50 

■32 . 4 7 

l . 90 

1 5 . 4 1 

10.00 

33.63 

7.11 

1 4 . 9 s' 

9,50 

3 4 . 70 

7 . 29 

15.1  3 

v.oo 

35.85 

7 . 36 

15. 7,s 

0 . 50 

3 7.00 

7 . 30 

15. 40 

ft , 0 0 

38  21 

7 . 1 2 

15.55 

*.50 

JV.  4 7 

1 . 90 

15.  0 

62 


APPENDIX  B 
Edited  Isotherms 

This  appendix  contains  the  entire  set  of  edited  time  series  of 
isotherm  depths  used  in  the  spectral  analysis  and  other  processing.  As 
mentioned  in  the  section  on  analysis,  only  series  that  were  at  least  80% 
complete  before  editing  were  chosen  to  be  edited  for  further  processing. 
Because  of  this,  not  all  tow  legs  are  represented. 
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APPENDIX  C 

Buoyancy  Frequency  Profiles 


87 


Plots  of  the  horizontally  averaged  buoyancy  frequency  are  presented 
in  this  section.  The  computational  method  is  described  in  Appendix  A. 
The  values  at  the  midpoints  between  thermistors  are  plotted.  All  plots 
use  the  same  scales,  which  are  shown  only  on  the  following  pane. 
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APPENDIX  D 

Spectral  Analysis  Methods 


Standard  spectral  analysis  techniques  were  used  in  processing  the 
data  for  this  study  and  only  pertinent  details  will  be  noted  here.  Ex- 
cellent references  for  the  subject  are  the  test  by  Bath  (1974)  and  the 
Ph.D.  thesis  by  Andreas+  (1977). 

The  complex  form  of  the  Fourier  coefficients  were  used  in  the  anal 

ysfs.  That  is,  a data  series  xk.  k - 0.  1 N-l  is  represented  as 

m_i  2nink 

1 " 1 ~1T~ 

Xk  " N kfQ  Xn  e for  k = 0,  1,  ....  N-l 

The  period  of  the  data  series  is  denoted  by  T and  the  sampling  interval 
by  A,  so  that  T = Na.  The  complex  Fourier  coefficients  are  given  by 


X 


n 


N-l  -2tTink 

nf0  xk  e N for  n = 0,  1,  ....  N-l 


Note  that  in  this  representation  *n  is  symmetric  about  f - 1,  if  X|(  is  a 
real  data  series.  That  is  X„_n  = Xn,  where  Xn  denotes  the  complex  con- 


jugate  of  Xp.  The  power  spectrum  is  then  given  by 


P 

xn 


? * 
= XX 

N2  n " 


for  n = 0,  1,  ...|  - i. 


The  power  spectral  density  is  defined  as  % - T Pxn.  if  Xfc  1s  a rea1 
series,  then  Pxp  are  also  real.  In  practice,  the  terms  "spectrum" 
and  "spectral  density"  are  often  used  interchangeaoly.  The  approach  taken 


Andrf?*,’J;.!:-,J977:  Observations  of  velocity  and  temperature  and 
e^tjmates  of  momentum  and  heat  fluxes  in  the  intpmai  Kn.m^a^v" 
laye£_over  Arctic  leads,  Ph.D. Thesis,  0re~gon  State  University 


in  this  study  is  to  use  "spectral  density"  but  to  refer  to  it  as  the 
"spectrum". 


Now  suppose  a second  data  series,  y^,  k = 0,  1 . ...  N-l  exists  with 
Fourier  coefficients  Yk-  The  cross  spectrum  of  the  two  series  is  given 
by 


7 * 

p = ^5-  X Y 
,2  n n 


for  n = 0,  1 , . . .^j-  - 1 


xyn  Nc 

and  cross  spectral  density  is  6 - T P . This  series  will  be  complex, 

A Jr  II  AJf  II 

with  the  real  part  called  the  cospectrum  and  the  imaginary  part  the  quad- 
spectrum.  The  equivalent  polar  representation  is  the  coherence  squared 
and  phase,  given  by 


Con 


(Co2  + Quad2) 


xyn 


rxn  Tyn 
Quad 

v - <^;n 

In  this  study  coherence  values  rather  than  coherence  squared  are  used. 

The  spectral  estimates  are  smoothed  by  bands  that  are  equally  spaced 
on  a logarithmic  scale.  Ensemble  averaging  of  several  series  is  done  in 
a similar  manner  of  the  coefficient  of  the  individual  series.  Confidence 
limits  can  be  calculated  by  assuming  a chi-squared  distribution  of  the 
estimates  in  each  band.  The  variance  is  used  to  compute  the  equivalent 
degrees  of  freedom  for  each  band  averaged  <j>:  v = varfance(,j,y-  The  true 
value  for  $,  denoted  as  <J>,  can  then  be  estimated  as  being  in  the  interval 


Xv:l-a 


2 

x 


v:a 


At  the  95%  confidence  level,  a = 0.025,  the  x values  can  be  approximated 
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4:0.025  “ * 9 892V-3'4)  (1  - |v  - l.OS&jV 

xv:0.975  ’ " fv  * '-96(|v'  * 

For  coherence,  a 95%  significance  level  is  computed  for  the  band  averaged 
values.  This  level  can  be  determined  after  Julian  (1975): 
y = (1  - t where  1-p  is  the  significance  level  (here  p is 

0.05  for  95%)  and  n is  the  number  of  samples  averaged. 

Various  schemes  for  windowing,  detrending,  and  filtering  the  series 
were  tried  and  abandoned.  Frankignoul  (1974)  studied  the  problems  of 
determining  internal  wave  spectra  from  short  series  and  found  end  point 
jump  to  be  a major  source  of  contamination.  One  method  he  suggested  to 
minimize  these  effects  was  chosen  for  the  present  study.  The  data  series 
is  first  pre-whitened  by  taking  the  first  forward  difference: 

yk  _ *k+i  “ *k  k - 0,  1 . . . , N-2 

yN-l  = yN-2 

After  the  spectrum  is  computed,  it  is  recolored  by  division  by 
2tt  k 

2(1  - cos  (^p)),  the  transform  of  the  differencing  scheme. 


1 
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APPENDIX  E 
Spectra 

Plots  of  the  autospectra  for  the  lowest  and  highest  isotherms  of 
individual  tow  legs  are  presented  in  this  appendix.  Again,  only  edited 
isotherm  series  that  were  at  least  80%  complete  before  editing  were  used. 
All  the  plots  have  the  same  scales,  which  are  given  on  the  following  page 
only.  Symbols  used  on  the  plots  are:  X for  the  spectrum  of  the  lowest 
isotherm  of  the  leg,  and  + for  the  spectrum  of  the  highest.  For  actual 
isotherm  values  and  depths,  see  the  listings  in  Appendix  A. 


RUN  1 3R  RUN  134- 
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RUN  1 7R  RUN  178 


(SW 


RUN  22K  RUN  22L 


